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Abstract 

Using (9.32,5.88) million T{2S,3S) decays taken with the CLEO III detector, we obtain five 
product branching fractions for the exclusive processes T(25) 7X60,i,2(l-P) — ^ ll'^i'^S) and 
T(35) — 7> 7X61,2 — ^ u'^i^S). We observe the transition Xboi^P) I'^i^S) for the first time. 
Using the known branching fractions for B[T{2S) 7Xbj(l-P)], we extract values for B[xbji^P) — ^ 
7T(1S')] for J = 0, 1, 2. In turn, these values can be used to unfold the T(35) product branching 
fractions to obtain values for B[T{3S) — > 7Xbi,2(l-P)] for the first time individually. Comparison of 
these with each other and with the branching fraction B[T{3S) "tXbo] previously measured by 
CLEO provides tests of relativistic corrections to electric dipole matrix elements. 



* Now at: Pacific Northwest National Laboratory, Richland, WA 99352 
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I. INTRODUCTION 



The bottomonium (bb) resonances display a rich pattern of electromagnetic transitions, 
including electric dipole (El) transitions between S'-wave (T{nS)) and P-wave {Xb{nP)) 
states Branching fractions for these transitions involving the lowest Xbi^P) states are 
summarized in Table | 0. Notable in Table | is the suppression of the transitions from the 
T{3S) to the X6j(lP) states. The electric dipole matrix element for the T{3S) — )■ 7X;,j(lP) 
transition, {lP\r\3S), is very small (see Ref. for a discussion), and thus is quite sensitive 
to assumed shapes of wave functions due to various relativistic corrections. In the non- 
relativistic limit it should be independent of J, with the corresponding decay rates given 
by 

T[T{3S) ^ 7Xw(lP)] = ^ « (2J + 1) E'^ \{lP\r\3S)\' . 

Thus, one would expect the rates for J = : 1 : 2 to be governed by the term x (2 J+1), or 
to be in the ratio 1:2.4:3.6. Various treatments of relativistic corrections to these decay rates 
1^ imply ratios differing considerably from these values and from one another, providing an 
opportunity to distinguish between them. 

TABLE L Previous data |||] for branching fractions for electric dipole transitions in the bottomo- 
nium (bb) system involving the lowest P-wave spin-triplet states Xbi^P)- Note the absence of 
measurements for T(35) — )• 7Xfei,2(l-P) and Xfeo(l-P) ~^ I'^i^S). 



Transition 


(MeV) 


B (%) 


T{3S) ^ 


■ 7X6o(l-P) 


483.9 


0.30 ±0.11 


T(35) - 


>7Xw(lP) 


452.1 


< 0.17 


T(35) ^ 


■ 7Xfe2(l-P) 


433.5 


< 1.9 


T{2S) ^ 


■ 7Xf)o(l-P) 


162.5 


3.8 ±0.4 


T(25) - 


>lXbi(.lP) 


129.6 


6.9 ±0.4 


T{2S) ^ 


■ 7Xfe2(l-P) 


110.4 


7.15 ±0.35 


XboilP) - 


7T(15) 


391.1 


< 6 


XbiilP) 


^ 7T(15) 


423.0 


35 ±8 


Xb2{lP) ■ 


7T(15) 


441.6 


22 ±4 



As seen in Table |, the four photon energies for the cascade transitions T(3S') — ?■ 
7Xfei,2(l-P) — 77T(1S') are all in the range 423-452 MeV in the rest frames of the de- 
caying particle, making it difficult to extract the product branching fractions for individual 
values of J Thus, until now, only the sum over J values (assumed here to be J = 1, 2) 

has been measured. For the summed product branching ratio, 

Bsura = Yl ^[^(3-^) ^ ^XbAlP)] X BixbAlP) 7T(1^)], 

J=l,2 

Ref. obtains -Bsum = (l-2^o 3 ± 0.09) x lO'^, while Ref. § obtains B^nm = (2.14 ± 0.22 ± 
0.21) X 10-3. 
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In this article, we obtain separate product branching fractions for J = 1 and J = 2. 
Notice, however, that unfolding B[T{3S) — j- 7x^1,2] from the product branching fractions 
requires knowledge of the X6i,2(l-P) — ^ I'^i^S) rates, which, as seen in Table |, have relative 
errors exceeding 20%. Therefore we also measure similar product branching fractions for 
T{2S) transitions, B[T{2S) 7X6j(lP)]x B[xbj{lP) ^ 7T(15)], from which the El 
rates for Xbj{^P) decays can be extracted by using the known T(2S') branching fractions 
in Table |. Thus, with product branching fractions for two-photon cascade transitions from 
T{3S) and T(25') to T(15') through Xbj(lP), we can make the first determinations of El 
branching fractions from T(35') through Xfej(l-P) for J = 1 and J = 2 and for Xboi^P) El 
transitions to T(15'), as well as improved values for Xbi,2(l-P) — ^ 7'^(1'S'). The branching 
fractions B[T{3S) — )■ 7x^1,2] may be compared with those for J = previously measured 
by CLEO [I, B[T{3S) 7X60(1^^)] = (0.30 ± 0.04 ± 0.10)%, providing tests of relativistic 
corrections to electric dipole matrix elements such as those involved in the predictions of 
Refs. i. 

We discuss the data samples in Sec. |T|. We describe our analysis method in Sec. |T|, 
and our fits to the transitions T(2S') — )■ 7Xbj(l-P) — )■ 77T(1S') and T(3S') — )• 7Xbj(l-P) 
77T(15') in Sees. |^ and ^ respectively. We summarize our results and compare theoretical 
predictions against the measured branching fractions B[T{3S) — )■ 7Xfej(l-P)] in Sec. and 
against B[xbj — > 7^^(15')] in Appendix]^. 



II. DATA SAMPLE 

Our event selection criteria (and this entire analysis procedure) closely follow the analysis 
of T(2S') riT{lS) The data used in this analysis were collected in e~^e~ collisions at 
the Cornell Electron Storage Ring (CESR), at center-of-mass energies at and ~ 25 MeV 
below the T(IS'), T(2S'), and T(3S') resonances. Integrated luminosities as well as the 
estimated numbers of these narrow resonance decays are shown in Table |T[ Events were 
recorded in the CLEO III detector, equipped with an electromagnetic calorimeter consisting 
of 7784 thallium-doped cesium iodide (Csl) crystals and covering 93% of solid angle, initially 
installed in the CLEO II [|l^ detector configuration. The energy resolution of the crystal 
calorimeter is 5% (2.2%) for 0.1 GeV (1 GeV) photons. The CLEO III tracking system ^ 



consists of a silicon strip detector and a large drift chamber, achieving a charged particle 
momentum resolution of 0.35% (1%) at 1 GeV/c (5 GeV/c) in a 1.5 T axial magnetic field. 



TABLE II. Integrated luminosities of data sets used in this analysis, in units of pb~ . "ON" 
corresponds to data sets taken in the vicinity of the nominal masses |^ of the corresponding 
narrow resonances, while "OFF" represents data taken ~ 25 MeV below the respective resonance 
masses. The estimated numbers of narrow resonance decays in each data set are listed at the 
bottom row. 



T{1S) T(25) T(35) 

ON 1056 1305 1387 

OFF 190 438 158 

#T(nS') (10^) 20.81 ± 0.37 9.32 ± 0.14 5.88 ± 0.10 
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III. ANALYSIS METHOD 



One can choose a set of ^^^"77 events = or fi^) in which the two photon energies 
sum up to a range consistent with the transition T(2S', 35*) 7Xbj(l-P) — 77T(1S'). We 
label the energy of the lower-energy photon E!^™, and that of the higher-energy photon Ejj's*^. 
In this section, we describe how we select our lepton candidates from T(IS') — j- £^i~ , how 
we apply kinematic constraints, and what our main backgrounds are. 



A. Selection of leptons and photons 

In order to identify leptonic decays of T(IS'), we first select the two highest-momentum 
tracks in an event. We call a track an electron candidate if > 0.75 or a muon candidate if 
E/p < 0.20, where E is the energy observed in the calorimeter shower associated with it and 
p is its momentum measured in the tracking system. Each track must satisfy | cos 6*1 < 0.83, 
where 6 is the angle with respect to the positron beam direction, to ensure reliable triggering 
and optimal performances of the tracking system and calorimeter for lepton identification. 
These tracks must originate within 5 cm (5 mm) along the beam direction (in the r-0 plane) 
of the the interaction point (IP). Both tracks must be of the same lepton type and be of 
opposite charge. 

Electron candidate tracks are dealt with somewhat differently, as they may radiate en- 
ergy via bremsstrahlung, and also contain significant contamination from radiative Bhabha 
scattering. To recover bremsstrahlung, we add to each lepton candidate's four-momentum 
the four-momentum of any photon candidates found to lie within a cone of 100 mrad of the 
lepton candidate track direction at the IP. To suppress contributions from Bhabha scatter- 
ing in 77e+e~ final states, we require e+ candidates to satisfy cos6'e+ < 0.5, where the final 
state positron makes an angle 6^+ with the incoming positron beam direction. This selection 
criterion greatly suppresses Bhabha scattering background while keeping a large fraction of 
the signal. 

Photon candidates must be detected either in the barrel (| cos6'| < 0.81) or in the endcaps 
(0.85 < I cos^l < 0.93) of the calorimeter. Each must have a lateral shower profile consistent 
with that of a photon, and the shower energy must exceed 30 (50) MeV in the barrel 
(endcaps). Additionally, such showers must not be aligned with the initial momentum of a 
track. 



B. Background composition 

As in the analysis of T(2S', 35") — 7T^/riT{lS,2S) 0, the dominant (and almost sole) 
sources of background are the doubly-radiative QED processes e'^e~ — )■ £+£"77. Such events 
can completely satisfy the restrictions on kinematic fit quality (see below) when the 
coincidentally has an invariant mass near that of the T(15'). Using our off-resonance data 
samples described in Table H to study these backgrounds, we find that such events produce 
smooth, nearly fiat, nonpeaking spectra in El^^ and Ejj^sh_ p^j. ^j-^g T(3S') — )■ analysis, 
we prepare Monte Carlo (MC) simulations of these processes for use in the nominal fits (see 

Sec. frx% ). 



To improve photon resolutions, we use values of kinematically constrained four-vectors 
instead of the observed ones. We perform kinematic fitting on events in the following manner 
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(cf. Ref. PI). We constrain the two leptons to have the T(IS') mass and the total four- 
momentum of the T(IS') candidate and the two photons to be that of the incoming e^e~ , 
accounting for the non-zero crossing angle of the beams. The dilepton vertex is constrained 
to the beam spot, which itself is measured with multi-track hadronic events for each run of 
approximately 1 hour duration. The photons are assumed to originate at the beam spot. 

The specific procedure employed is the following: Obtain reduced values (i.e., values 
of per deg ree of freedom) from the above vertex and momentum fits. Call them 'x^ 
and Xm ic; respectively. Combine the mass-fitted object with two-photon candidates and 
constrain the sum of their four-momenta (a 4C fit) to the laboratory four-momentum, ob- 
taining a reduced value denoted by Xm4C- Require reduced values in the following 
order: (1) xlw < 10; (2) xl,w < 10 (3) xl,,c < 5- 

Our kinematic fitting software assumes that measured distributions are Gaussian, but 
there are well-understood and reasonably well-modeled low-side tails on the energy response 
of the calorimeter. Hence, the (reduced) distributions from fits in both MC and data 
will also have tails not seen in formal probability distributions. Good agreement of the 
reduced x^ distributions between data and MC is essential. We explored this agreement 
in previous CLEO analyses [jl2| and in this analysis (Sec. |1 V C 4|) . We included the small 
discrepancies we found in our systematic errors. 

Figure |l]shows distributions of one of our reduced x^ variables, Xm 4C' ^'^e T{2S, 3S) — t- 
77T(1S') analyses. Standard and alternate restrictions on these variables are indicated in the 
figure. The overlaid histograms indicate the contributions of MC signal and background, 
weighted by our final measured values. Based on the scaled off-resonance data, the ex- 
pected background levels in the signal regions correspond to background-to-signal ratios of 
~ 1.5 (3)% for the T{1S) candidates and ~ 15 (20)% for the T{1S) e+e" 

candidates in T(25') (T(35')) decays. 



IV. ANALYSIS OF T(25) -y-fl+i- VIA Xbj{lP) 

Because the expected E!^"^ spectra (110-160 MeV) have excellent separation of decays 
through different Xfej(l-P) spin states, we can obtain the yields of the individual J contribu- 
tions with a fit to Ely°™, without regard to Elj^sh (390-440 MeV). The El^°" spectra exhibit three 
clearly distinguishable peaks with known peak energies and resolutions entirely dominated 
by measurement effects, although detector resolution is improved upon by the constrained 
fit. Figure |2| illustrates the E?^" distributions from data, along with the fits described in the 
following subsections. 



A. Signal Monte Carlo samples 

We use the EVTGEN event generator [T^ to generate expected signal shapes in Monte 
Carlo (MC) simulations. To simulate the photon angular distributions, we assume that the 
radiative transitions are pure electric dipole. We generate 100 k events with ;B[T(25') — )■ 
7Xw(lP)] = B[xbj -> 7T(1^)] = 100% while B[T{1S) ^ e+e"] = B[T{1S) -> fi+fi-] = 
50% for each of J = 0, 1, and 2. We fit each of the MC distributions of E^°^ with a double 
Gaussian whose difference between the two means we allow to float. Resultant reconstruction 
efficiencies are shown in Table |IT|. 
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2541110-008 




|2 



5Cm,4C 

FIG. 1. Distributions of Xm4C ^"^^ T(IS') /^'''/i" candidates (a,b) and T(15') e+e" candidates 
(c,d) based on events of T{2S) jjT{lS) candidates (a,c) and T{3S) jjT{lS) candidates 
(b,d) via Xbji^P), respectively. Open circles represent data, solid histograms the sum of scaled 
signal MC and background contributions, where the background levels are indicated by scaled 
off-resonance data (dotted histogram) and, for T{3S) decays, QED MC simulation (dashed his- 
tograms). Solid arrows indicate standard selection criteria, and dashed arrows alternate values 
used for systematic studies of dependence on selection criteria. 



B. Fitting the data 

Our nominal fit procedure is to take these double Gaussian shapes based on the signal 
MC samples to fit to data, fixing the respective narrower Gaussian widths and differences 
between the two Gaussian means but allowing the larger widths {aj where J = 0, 1, and 2) to 
float. We use a flat background shape whose normalization is also allowed to float. We then 
perform a maximum likelihood fit. Figure |] shows fits to data for T(15') — j- candidates 
(a,b) and T(IS') -> e~^e~ candidates (c,d) with flat background shapes represented by the 
dashed histograms, while the dotted histograms are based on scaled off-T(2S')-resonance 
data. In Fig. |2|(b), we zoom in to a smaller vertical scale to emphasize the J = component, 
which is clearly visible for the fj,^fj.~ candidates. However, in Fig. ||(d) the J = peak is 
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2541110-007 




70 90 110 130 150 170 190 



e;™ (MeV) 

FIG. 2. Fits to data with a flat background shape represented by dashed histograms for T(15') — >■ 
fj.'^fj-" candidates in the two plots (a,b) [(c,d)]. Plots (b) and (d) are identical, respectively, 

to (a) and (c) but zoomed in so as to highlight the J = component. The dotted histograms 
represent the scaled background shapes based on the off'-resonance data. The dot-dashed lines 
show components of each of the three photon lines. 

obscured by the relatively larger background in the e+e~ candidates. To compensate for 
the larger backgrounds in the distributions for the the e~^e~ candidates, we fix the ratio 
of the widths of the J = peak to the width of the J = 1 peak (o"o/cri) to be equal to 
the ratio obtained for the signal MC samples. (In our systematic study, we remove this 
restriction and observe the deviation from our central value.) The observed yields, along 
with efficiency-corrected products of branching fractions, are shown in Table |ITT[ The last 
column of Table |IT| (labeled as "All J" ) represents measurements of these same two photon 
cascade events, but is summed over J for J = 0, 1, and 2. We obtain the yield simply by 
subtracting the fitted background shape, with normalization determined from the nominal 
fit, and then by summing the resultant spectrum over the signal region. There, efficiencies 



8 



TABLE III. Fitted yields, reconstruction efficiencies (e), and corresponding branching fractions are 
shown with statistical errors only. Here, Bl = B[T{2S) 7Xbj(lP)], B2 = B[xbjilP) ^ 7^(15')], 
and B3 = B[T{1S) £+£"]. The last column ("Ah J") shows resuhs for the sum of J = 0, 1, and 
2 components obtained by subtracting fits to the backgrounds from the data, as described in the 
text. 







J = 


J= 1 


J = 2 


Ah J 


Yields 




16 ±9 


1068 ± 36 


600 ± 30 


1684 ± 41 


Yields 




71 ± 10 


2154 ± 50 


1170 ± 39 


3395 ± 58 


e (10-2) 




20.9 ±0.2 


21.3 ±0.2 


20.2 ±0.2 


20.9 ±0.2 


6 (10-2) 




38.6 ±0.3 


39.9 ±0.3 


37.7 ±0.3 


39.1 ±0.2 


BlxB2x B3 (10-4) 




0.083 ± 0.044 


5.38 ±0.18 


3.19 ±0.16 


8.66 ±0.21 


BlxB2x B3 (10-4) 




0.196 ±0.028 


5.79 ±0.13 


3.33 ±0.11 


9.32 ±0.16 


BlxB2x B3 (10-4) 


e+e- and 


0.163 ±0.024 


5.65 ±0.11 


3.29 ± 0.09 


9.08 ±0.13 



are weighted by the measured branching fractions for each spin state. The values for the 
fits are 50.1 for 60 data points (minus 9 parameters), c.L = 51.0% for T(15') — > e~^e~ , and 
51.4 for 60 data points (minus 10 parameters), c.L = 42.0% for T(IS') — ■ 



C. Systematic uncertainties 



Systematic uncertainties are assessed as J-independent contributions that impact the 
yields of several modes equally, or are determined individually. Where errors differ for e+e" 
and /i"*"/!- candidates, they are averaged with the same statistical weights used for combining 
the respective product branching fractions. 



1. Uncertainties common to both e"*"e and ^'^ jjL candidates 

The uncertainties common to both e^e~ and fi^fi" candidates are relative uncertainties 
known already from external sources, and which impact every yield identically. The relative 
uncertainties on the numbers of resonance decays in our data sets were estimated in Ref. 



p. We take the dilepton reconstruction systematic uncertainty from Ref. [|14| as 1.0%, since 
our dilepton reconstruction is identical to that used there. Similarly we use the result from 
that paper for vr'^-finding, which was 1.6% per 7r°. We therefore use this as an estimate of 
the uncertainty for reconstructing the two photons in the cascade, since the photon energies 
we are studying in this analysis very closely resemble the photon energies in Ref. [|14|, and 



otherwise the processes are kinematically very similar. Instead of 1.6% for the pair, we 
round up to 2.0% to be conservative. 
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2. Signal shapes 

In our nominal fit procedure, we take a double Gaussian, fitted to the signal MC samples, 
to represent the signal line shapes but float the larger widths of the Gaussians for the fits 
to data to accommodate imperfect simulations of detector resolutions. Here, we have tried 
the following variations: 

• Constrain all aj in both T(IS') — £'^£~ candidates while requiring the ratios ao/ai 
and o^jox to be the same as in our signal MC samples. 

• Do not constrain the "a" in either lepton flavor. 

3. Background shape 

The only backgrounds predicted by a MC simulation of all T{2S) decays are from 
T{2S) — 7- 7r°7r°T(lS'). Even so, this background source is found to be negligible because 
it contributes only 5-10 events to the signal regions and it has no significant structure in 
photon energy. Hence, to the extent it matters at all, it will tend to get absorbed into the 
background shape in the fit. Backgrounds from T(2S') — )■ 7i^/riT{lS), using the recently 
measured branching fractions |^, are also found to be insignificant. 

To represent our background, whose main compositions are either doubly-radiative 
Bhabha events or //-pairs, we use a fiat shape with floating normalization in our nom- 
inal fit. To probe the sensitivity to the fitted yields due to this assumption, we try a 
first-order polynomial. 

To determine possible systematical effects due to the fixed background normalization 
in the "All J" procedure, we vary the background by ±1 cr of the statistical uncertainty 
observed in our nominal procedure. 

4- Kinematic fit reduced requirements 

We consider different choices in reduced criteria for our vertex and four-momentum 
fits. The four variations to our standard selection are: XmAC < 4 or < 10 rather than 5, 
Xmic < 5 rather than 10, and xlic ^ rather than 10. The spin of the Xbj{^P) state 
should have almost no influence on how these reduced distributions are simulated, so to 
increase statistical stability on the resultant variations, we take variations from the "All J" 
procedure to assign this systematic uncertainty. 

5. Lepton flavor difference 

We have also assessed possible systematic uncertainties due to the difference between 
e~^e~ and ^~ results by calculating the yield for the sum of J = 1 and J = 2 for e^e~ and 
then for yU"*"/!". (As in the case of determining possible uncertainties due to requirements 
on the fitted reduced x^, the spin of the Xbji^P) state should have no influence on whether 
leptons are correctly reconstructed.) We then took half the difference between the yields 
obtained for each lepton flavor, and divided by the average yield as an estimate of the relative 
systematic uncertainty arising from lepton flavor differences. 
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6. Additional contributing uncertainties 



Other possible systematical effects we have investigated include: variations in fit ranges, 
histogram binning, statistical uncertainties in signal MC samples, uncertainties in the mea- 
sured B[T{2S) 7Xw(lP)]x B[xbj{lP) ^ 7T(15)] to weight efficiencies for the "All J" 
well as cuts on E/p to identify lepton species. 

Table |^ shows sources of systematic uncertainties we have considered for this analysis. 



TABLE IV. Fractional uncertainties (in %) on the combined-dilepton product branching fractions, 
B[T{2S) 7X6j(l^')]x l3[xbj{lP) 7T(15)]x B[T{1S) £+£-], due to variations of the listed 
selection criteria and fit procedures. Tfie last column ("All J") shows fit results for the sum of the 
three two-photon cascades via Xbj{^P) states for J = 0, 1, and 2. For this last column, the entry 
for "MC simulation" includes not only statistical errors on reconstruction efficiencies, but also the 
effect of total uncertainties of the measured branching fractions on the weighted efficiency. 



Contribution 


J = 


J = 1 




J = 2 


All J 


-^T(2S) 






1.5 






Track-finding 






1.0 






Photon-finding 






2.0 






Reduced requirement 






1.4 






Lepton identification 






0.4 






Lepton flavor difference 






3.6 






Fit range 


7.1 


0.2 




0.7 


0.4 


Signal shape 


1.4 


0.5 




0.6 




QED bkg shape 


1.9 


0.01 




0.2 


0.9 


Bin width 


1.6 


0.04 




0.01 


0.01 


MC simulation 


0.8 


0.8 




0.8 


0.6 


Total 


9.0 


4.8 




4.9 


4.9 



D. Results on analysis of T(25) -yyi+i- 

Table shows our final results for our T(2S) analysis as well as those from other exper- 
iments. 



V. ANALYSIS OF T(35) 77£+£- VIA Xbj{lP) 

In the case of the three transitions T(3S)— )-7X6j(l-P)— )-77T(lS), the higher and lower en- 
ergy photons have similar energies. Furthermore, among these transitions in some cases the 
higher energy photon is emitted from the T(3S), and in some cases the higher energy photon 
is emitted from the boosted Xbj{^P)- Therefore, instead of fitting just one distribution, as 
we do for the T(2S) analysis, we maximize our use of information and fit the two-dimensional 
histogram of Ej^™ vs. E!^^^*^, where Ej^™ < Ejj'^'^ are kinematically constrained E^. Our fit 
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TABLE V. Final results of this analysis. Here, Bl = B[T{2S) 7Xw(lP)], 132 = B[xbj{lP) ^ 
7T(15)], and B3 = ^[T(15) ^ We use B3 = (2.48 ± 0.05)% § and Bl values from Table | 

to extract Bl x B2 as well as B2. The last column ("All J") shows fit results for the sum over 
J = 0, 1, and 2. Again, the first errors are statistical, the second errors are systematic, and 
the third errors (when applicable) are uncertainties due to uncertainties in Bl and/or B3. In the 
bottom half of the table we also show results from other experiments for a comparison. 





.7 = 


J = 1 


J = 2 


All J 


ei X B2 X B3 (10-*) 
61 X B2 (10-=*) 
B2 (10-2) 


0.163 ±0.024 ±0.015 
0.659 ± 0.096 ± 0.059 ± 0.013 
1.73 ±0.25 ±0.16 ±0.19 


5.65 ±0.11 ±0.27 
22.8 ±0.4 ±1.1 ±0.5 
33.0 ±0.6 ±1.6 ±2.0 


3.29 ±0.09 ±0.16 
13.3 ±0.4 ±0.6 ±0.3 
18.5 ±0.5 ±0.9 ±1.0 


9.08 ±0.13 ±0.44 
36.7 ±0.6 ±1.8 ±0.7 



Values of B2 (10 ^) from other experiments 



PDG average [2] n < 6 at 90% CL 35 ± 8 22 ± 4 — 

Crystal Ball [15] LJ < 6 at 90% CL 32±6±7 27±6±6 — 

CUSB CH < 11 at 90% CL [16] 47 ± 18 [17] 20 ± 5 [17] — 

□ □ □ 



will utilize the 2D histograms from our signal Ivlonte Carlo for each of the J = 0, 1, and 2 
samples, as well as MC samples of doubly-radiative Bhabha events and /i-pairs. 



A. Monte Carlo samples 

1. Signal MC for T{3S)^-fXbj{lP) 

Approximately 100k events for each spin were generated for each of the Xbj(l-P) subsam- 
ples, where T(3S)^7Xbj(lP) at 100%, Xbj{lP)^l^{lS) at 100% and with T(1S) decaying 
half the time to each of fi'^fi~ and e+e~. These samples were used to generate the three 
relevant 2D histograms of E!^"™ vs. E^^^^ for each spins. In projections of these histograms 
on the E^^^^ and E!^"™ axes, for J = and J = 1, the higher-energy photon has a sharp 
distribution, while the lower-energy photon is Doppler-broadened (see Table |). For J = 2, 
the situation is reversed to a large extent. 

Figure ^ shows the Elj^™ vs. Eij'^^ two-dimensional histogram for T(35') signal MC samples, 
weighted according to our final measured branching fractions. Notice that the events are 
restricted to the diagonal band by two simple kinematic facts: first, that E)^^ < Elj^^*^, and 
second, that the invariant mass recoiling against the two photons has been constrained to 
that of the T{1S). The lower boundary line corresponds to the case in which the photon 
three-momenta have zero opening angle in the initial e~^e~ rest frame (recoiling against the 
dilepton pair with maximum possible combined momentum), and the upper boundary line 
to the case in which the photon momenta directly oppose one another (with the softer one 
traveling along the same direction as the T(IS') candidate, and the more energetic one in 
the opposite direction). Note that the boundaries are purely kinematic in nature due to the 
T(IS') mass constraint, and hold for background as well as signal events. In Fig |^, these 
limits are labeled 9.^.^ = and tt, respectively. 

The population within the band is directly related to the cosine of the angle 6^^ between 
the two photons in the rest frame of the X6j(lP). The distributions Wj{cos9^^) for J = 1, 2, 
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2541110-004 



435 



> 

CD 



-5 ^ 
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(MeV) 



FIG. 3. Distributions of E!^°* vs. e!^'^^ based on T(35) signal MC samples for J = (open circles), 
J = 1 (open rectangles), and J = 2 (closed triangles) using /i selection. The diagonal band, 
edged by lines labeled as Q-y-y = vr (or 0), where Q^^ is an opening angle between the two emitted 
photons, is generated due to our kinematic constraints. The three samples are normalized to our 
measured production rates in this work while the size of a symbol in a bin is proportional to the 
number of events for the corresponding J in that bin. 
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normalized so that their integral over cos is 1 , are 

15 1 

Wi{cose^^) = —(1 + -cos^^^^) , (1) 
iy2(cos^^^) = ^(73 + 21 cos' d^^) , (2) 

implying a slight enhancement at each end of the cos 9^^ range and hence at maximum and 
minimum Doppler-broadened photon energy. This feature is present in all of our signal 
MC samples that take account of the photon angular distributions properly. The respective 
reconstruction efficiencies are shown in Table |V^ . 

TABLE VI. Fitted yields, reconstruction efficiencies (e), and corresponding branching fractions are 
shown. Here, Bl = B[T{?>S) 7Xw(l-P)], B2 = B[xbj{lP) ^ 7T(15)], and B3 = B[T{1S) 

The last column ("J = 1 and 2") shows results of sum of the two two-photon cascades 
via Xbji^P) states for 1 and 2 obtained by subtracting fits to the backgrounds from the data, as 
described in the text. 





i+r 


J = 1 


J = 2 


J = 1 and 2 


Yields 




12 ±6 


48 ±9 


61 ±8 


Yields 




38 ±9 


78± 11 


117 ± 11 


e (10-2) 


e+e- 


22.3 ±0.2 


21.3 ±0.2 


21.6 ±0.2 


e (10-2) 




41.1 ±0.2 


38.8 ±0.2 


39.5 ±0.2 


BlxB2x B3 (10-5) 




0.91 ±0.49 


3.88 ± 0.70 


4.79 ±0.61 


BlxB2x B3 (10-5) 




1.58 ±0.38 


3.40 ± 0.49 


5.06 ± 0.47 


BlxB2x B3 (10-5) 


e^e~ and 


1.33 ±0.30 


3.56 ±0.40 


4.96 ± 0.37 



2. Background from QED processes 

Only the processes e'^e~^^^{e'^e~ , /x^/i-) (and to a much lesser extent, T(3S)— 7'7r°7r°T(lS)) 
can be reasonably expected to contribute much background. To represent the doubly- 
radiative QED events, we prepare MC samples generated via the Babayaga event genera- 
tor |18[ with sizes of roughly 100 (200) times larger luminosity than what the data have 
for doubly-radiative Bhabha (/i-pair) events, respectively. The very large number of QED 
events are then reduced with a very loose selection at the generator level to have at least 
two photons and a dilepton invariant mass near that of the T(IS'). Only those chosen events 
are passed along to the next stage of processing, the CLEO detector simulation. 



B. Fitting the data 

Figure ^ shows a distribution of El^°™ vs. El^J'^^ based on the on-T(35') data. To constrain 
the backgrounds in our fits, we use fit ranges for El^^ and Eij'^^ larger than the ranges 
illustrated in the figure. We choose our fit ranges to be (420 <E!J's'^< 560) MeV and 
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2541110-009 

450 I 1 1 1 1 1 1 1 1 1 1 1 1 




425 435 445 455 465 475 485 



Ej'^' (MeV) 



FIG. 4. Distributions of EJ^™ vs. E^'^ based on the on-T(3«S') data. Open triangles (closed 
rectangles) represent data points based on e~'"e~ selection. The size of each symbol is 

proportional to the number of events in the bin. 



(340 kE^^'^k 460) MeV. While the upper (lower) bound of E':^°"(E!j'g^) does not matter much 
because of bins containing no events, the choice of the lower (upper) bound of E1^^(E!J'^'^) 
controls the statistics available for fitting the backgrounds. Since the minimum El^^ is 
related to the maximum E^^^*^ by the kinematic constraints, we study the variation of the 
background scale factor as a function of the maximum E^^^^. Based on this exercise, we 
choose the maximum E!J'§^ to be 560 MeV and the minimum El;'^= 340(~ 900 - 560) MeV 
where E^^^*^ + El^" ~ 900 values for which the fitted normalization scale factors become 
stable compared to the sizes of their statistical errors. 

Using the QED MC background histograms as the background function and the three 
signal Monte Carlo samples, and fixing the normalization of the J = component with 
the measured B[T{3S) — )■ ■jXwi^P)] and B[xbo{^P) 7'^(15')] (this work), we perform 
a maximum likelihood fit to the the 2D data distribution in E!^°" vs. E^'sh_ flo^t the 
normalizations of both background and signals (for J = 1 and 2 only). The results of the 
fits projected onto the two photon-energy axes are shown in Fig. |^. 

We also extract signal yield by subtracting the fitted background shape, with normal- 
ization fixed based on the above nominal fit procedure, and then by summing the resultant 
distribution over the signal region. This yield corresponds to the one from both transitions, 
T(3S') — > 'j'j£~^£~ via Xbj{^P) for J = 1 and 2. The observed yields, along with efficiency- 



corrected branching fractions, are shown in Table VI 
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2541 1 1 0-006a 




380 400 420 440 435 455 475 495 



E{°^ (MeV) E,''^' (MeV) 

FIG. 5. Projections from the nominal fit to data onto (a,c) EJ^™ and (b,d) eJJ'^'^ axes for data 
using the (a,b) and (c,d) e+e~ selections. The dashed histograms represent signal photons 

via Xbji^P)- The dotted histograms represent the QED MC sample contribution as normalized by 
the fit, and the solid histogram is the sum of background and signal histograms. 



C. Systematic uncertainties 



We study the systematic uncertainties in a similar way as in Sec. IV C. We mention only 



studies whose methods differ from the description given in Sec. |IV Q . Table [VII| shows a 
summary of estimated systematic uncertainties. 



1. Binning scheme 

We calculate the branching fractions using 4 MeV bins instead of the nominal 2 MeV 
bins. Besides a difference in bin sizes, we also consider a difference in bin shapes with a 



16 



distribution of ^sum (= + E^^^^) vs. ^diff (= E^'^^ - E^^"^). 

2. Photon energy resolution 

In order to account for the possibility that the width of the signal MC photon distributions 
may not accurately reflect the width of the photon distributions in data, we have generated 
additional signal MC sets with the detector resolution broadened and narrowed by 10% of 
itself and recalculated the branching fractions using these altered MC sets. 

3. Photon absolute energy calibration 

To estimate possible systematic effects due to miscalibration of the absolute photon en- 
ergy, we examine how well the peak values in the kinematically-fitted distributions of El^^ 
for T{2S) — )■ 77T(1S') candidates are calibrated with respect to those expected. The largest 
deviation is 0.2 MeV for the J = 2 peak position when an e~^e~ final state is selected. It 
is well reproduced in our MC samples with input masses of T(1S', 2S) as well as of Xfe,/(1-P) 
based on the latest information 0. These photon energies are kinematically-fitted vari- 
ables, with EJ^" and Ei^'^h constrained to a fixed value for any given T(IS') momentum, so 
any shift in El^™ must be compensated by a corresponding one of opposite sign in E^^^^. 
Based on these observations, we conservatively vary signal positions by E!J'^^±1 MeV and 
simultaneously by E}!^^^1 MeV. 

4- Additional contributing uncertainties 

Other possible systematic effects that we have investigated include a variation in the 
(small) branching fraction for the transition through Xboi^P) and insertion of an explicit 
fixed background component from T(35') — )■ 7r°7r°T(lS'). 

D. Results on analysis of T(35) 77£+£~ 

Taking the above systematic uncertainties into account, we now arrive at the final results 
for the product of branching fractions for each of the transitions as shown in Table [Vlllj . 
Also shown are comparisons to results from other analyses. The first uncertainty is the 
statistical uncertainty, the second is the overall systematic uncertainty, and the third (when 
applicable) is the uncertainty due to external inputs. 

VI. RESULTS AND CONCLUSIONS 

We obtain product branching fractions for the exclusive processes T(2S') — )■ 7X6o,i,2(l-P) ~^ 
77T(15) (Table 0) and T(35) ^ iXbiM^) ^ ll'^i'^S) (Table IVTT]), where T{1S) is iden- 
tified by its decay to e+e" and ■ 

The extracted B[xbj{)-P) — ^ 7^(1*5')] are the most precise to date for J = 1, 2, while 
for J = this represents the first observation of this transition. These branching frac- 



17 



TABLE VII. Estimates of relative systematic uncertainties (in %) for this analysis, for each J, 
and the sum of J = 1 and 2. The first six entries are calculated for both spins. The fit to the 
sum of J = 1 and 2 has better stability against our binning scheme variations since yields of 
J = 1 and 2 are statistically anti-correlated, resulting in smaller variations in terms of their sum. 
For the rightmost column, the entry for "MC simulation" includes not only statistical errors of 
reconstruction efficiencies, but also the total uncertainties of the measured branching fractions on 
the weighted efficiency. 



Contribution 


J = 2 


J= 1 


J = 1 and 2 


-^T(35) 




1.7 




XI acK-iiiitiing 




1 n 




Photon-finding 




2.0 




Reduced requirement 




4.9 




Lepton identification 




1.1 




Lepton flavor difference 




1.2 




QED background 






8.8 


Include T(3S)^7r07rOT(lS) 


0.3 


1.2 


0.5 


Xfeo(l-P) Yield variation 


0.0 


0.6 


0.2 


Binning scheme 


3.3 


11.6 


2.7 


Photon energy resolution 


2.3 


3.7 




Photon absolute energy 


9.1 


11.0 




MC simulation 






0.7 


Total 


11.6 


17.5 


11.0 



TABLE VIII. Final results of this analysis. Here, Bl = B[T{3S) 7Xbj(lP)], B2 = B[xbj{lP) ^ 
7T(15)], and B3 = B[T{IS) £+i~]. We use B3 = (2.48 ± 0.05)% g and B2 = (18.5 ± 0.5)%), 
(33.0 ± 0.6)% from Sec. for J = 2 and J = 1, respectively, to extract Bl x B2 as well as Bl. 
The first three rows show results from this work while the second three rows show previous results. 
Here, the first errors are statistical, the second errors are systematic, and the third errors (when 
applicable) are uncertainties due to external sources. 







J = 1 


J = 2 


J = 1 and 2 


SI X S2 X S3 (10-5 


) 


1.33 ±0.30 ±0.23 


3.56 ±0.40 ±0.41 


4.96 ±0.37 ±0.55 


SI X S2 (10-4) 




5.38 ± 1.20 ±0.94 ±0.11 


14.35 ±1.62 ±1.66 ±0.29 


19.99 ± 1.50 ± 2.20 ± 0.40 


SI (10-3) 




1.63 ±0.36 ±0.28 ±0.09 


7.74 ± 0.88 ± 0.88 ± 0.38 




SI X S2 X S3 (10-5) 


[6] 


D 




5.20 ±0.54 ±0.52 


SI X S2 (10-*) [19] 
SI (10-3) [7] ^ 


□ 


< 1.9 


< 20.3 (11±6±2±1) 


I2tt ± 0.9 
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TABLE IX. Comparison of measurements and theoretical predictions for suppressed El tran- 
sition rates r[T(35) ^ 7Xw(lP)] and ratios Tj=i/Tj=o = r[T(35) ^ 7Xw(li^)]/r[T(35) ^ 
-fXboilP)], rj=2/rj=o = r[T(35) ^ 7X62(l^)]/r[T(35) ^ IXboilP)], and rj=2/rj=i = 
r[T(35) ^ 7Xb2(l^)]/r[T(35) ^ 7X6i(lP)]. The CLEO III values are based on rtotai[T(35)] = 
(20.32 lb 1.85) keV and are obtained by taking the central value of the measurement for the 
J = state ^] and the values for J = 1 and 2 from this work. The last row shows Tj=i/Tj=q and 
rj=2/rj=o when scaling rates according to x (2J + 1). 





rj=o (ev) r 


j=i (eV) rj=i/rj=o 


rj=2 (eV) Tj-_ 


=2/rj=o 


rj=2/rj=: 


CLEO III (This expt.) 




33 ±10 


0.54 ±0.25 


157 ±30 2.5J 


B± 1.01 4.75± 1.7f 


Inclusive expt. Q 


61 ±23 












Xbji^P) exclusive expt. Q 


< 186 


< 38 




< 413 






Moxhay-Rosner (1983) 


25 


25 


1.0 


150 


6.0 


6.0 


Gupta et al. (1984) 


1.2 


3.1 


2.6 


4.6 


3.8 


1.5 


Crotch et al. (1984) (a) 


114 


3.4 


0.03 


194 


1.7 


57 


Crotch et al. (1984) (b) 


130 


0.3 


0.002 


430 


3.3 


1433 


Daghighian-Silverman (1987) 


42 


(c) 


(c) 


130 


3.1 


(c) 


Fulcher (1990) 


10 


20 


2.0 


30 


3.0 


1.5 


Lahde (2003) 


150 


110 


0.7 


40 


0.3 


0.4 


Ebert et al. (2003) 


27 


67 


2.5 


97 


3.6 


1.4 


X (2 J + 1) 






2.4 




3.6 


1.5 



(a) Scalar confining potential, (b) Vector confining potential. 



(c) The authors did not provide a prediction for r[T(3S') ^Xbii^P)]- 



tions appear to be systematically smaller than the theoretical predictions (see Appendix ^), 
indicating that the hadronic widths of Xbj{^P) might have been underestimated. 

The extracted ;B[T(35') — ?■ 7Xfei,2(l-P)] may be compared with the branching fraction 
previously measured by CLEO, i3[T(3^) ^ 7Xw(l^)] = (0.30 ±0.04 ±0.10)% 0, providing 
tests of relativistic corrections to electric dipole matrix elements. Table shows comparison 
against some theoretical predictions in terms of transition rates as well as ratios of transition 
rates while Fig. ^ shows the ratios pictorially. It might be worth revisiting these calculations 
in light of our new experimental results. 
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FIG. 6. Illustration of ratios of suppressed El transition rates, Tj=i/rj=Q (circles), rj=2/^j=o 



(squares), and rj=2/^j=i (triangles) from Table IX. 



Appendix A: Comparison of branching fractions B[xbj{^P) 7'^(1'S^)] with pre- 
dictions 

The measured branching fractions for Xbji^P) — >" I'^i^S) may be compared with the 
predictions ^ summarized in Table ^ Most of the predicted branching fractions for these 
electric dipole transitions are systematically larger than the experimental values, indicating 
that the hadronic widths Th were underestimated. A modest increase in the assumed value of 
as{ml) leads to much better agreement with experiment. As one example, the values in Ref. 
Il20| were calculated for as{rnl) = 0.18. For this value it was found that Th[xb{o,i,2){^P)] = 
(791,38.3,132.3) keV, while the El transition rates were predicted to be T[xb(o,i,2){^P) 
7T(1S')] = (26.1,32.8,37.8) keV. The hadronic widths scale for the J = and J = 2 states 



ag{mi) times known QCD correction factors while for the J = 1 state they scale as 



Q;|(m^). (The QCD correction factor for J = 1 is not known pT| and will be ignored.) 

Using the scale factors, values of a5, and the QCD correction factors described in the 
previous paragraph, we predict 

ThlxmilP)] _ (cis{ml)V 1 + lQ.Qas{ml)/T, 



791 keV V 0.18 / 1.573 

F,[xw(lP)] _ (as{ml) , 



38.3 keV V 0.18 

F,[x.2(lP)] _ ( o^s{ml) V l-Q.las{ml)/n 

132.3 keV V 0.18 ) 0.994 ' ^ ' 

The minimum x^ of a fit to the values determined in this work is found to be 1.45 for 2 
degrees of freedom The value of as that minimizes x^ and the corresponding error (defined 
by the range for which Ay^ < 1 from the minimum) are as'(m^) = 0.214 ± 0.006. This is 
quite consistent with the determination of Ref. |^ for a scale of 5 GeV (see Fig. 5 there). 



20 



TABLE X. Comparison of our results for B[xbji^P) 7"f (15')] with some theoretical predictions 
|, 1^, in units of 10"^. 



Reference 


J — u 


T 1 

J — 1 


7 O 

J — / 


CLEO III 


1.73 ±0.35 


33.0 ±2.6 


18.3 ± 1.4 


Moxhay-Rosner (1983) 


3.8 


50.6 


22.3 


Gupta et al. (1984) 


4.1 


56.8 


26.7 


Grotch et al. (1984) (a) 


3.1 


41.9 


19.4 


Grotch et al. (1984) (b) 


3.3 


43.9 


20.3 


Daghighian-Silverman (1987) 


2.3 


31.6 


16.6 


Kwong-Rosner (1988) 


3.2 


46.1 


22.2 


Fulcher (1990) 


3.1 


39.9 


18.6 


Lahde (2003) 


3.3 


45.7 


21.1 


Ebert et al. (2003) 


3.7 


51.5 


23.6 



(a) Scalar confining potential, (b) Vector confining potential. 



At this value, the rescaled values predicted in the approach of Ref. |20| are B[xbo,i,2i^P) 
7T(1^)] = (2.1, 33.8, 16.8)% and Ttot[Xbo,iMP)] = (1221, 97, 225) keV. 
The predicted ratio 

r,(x.o(lP)) 15 1±10.0«5(mg)/7r 
- r,(x62(lP)) 4 l-O.lasimD/n ^ > 

is roughly R = 5.92 ± 12[a5(m^) — 0.18], to be compared with the value R = 8.6 ±3.2 based 
on the observed branching fractions. Thus, the QCD corrections go in the right direction to 
modify the uncorrected value of 15/4 = 3.75. 
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